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Abstract 

In  our  prior  work,  we  presented  a  highly  effective  lo¬ 
cal  search  based  heuristic  algorithm  called  the  Largest 
Expanding  Sweep  Search  (LESS)  to  solve  the  minimum 
energy  broadcast  (MEB)  problem  over  wireless  ad  hoc 
or  sensor  networks.  In  this  paper,  the  performance  is 
further  strengthened  by  using  iterated  local  optimiza¬ 
tion  (ILO)  techniques  at  the  cost  of  additional  compu¬ 
tational  complexity.  To  the  best  of  our  knowledge,  this 
implementation  constitutes  currently  the  best  perform¬ 
ing  algorithm  among  the  known  heuristics  for  MEB.  We 
support  this  claim  through  extensive  simulation  study, 
comparing  with  globally  optimal  solutions  obtained  by 
an  integer  programming  (IP)  solver.  For  small  net¬ 
work  size  up  to  20  nodes,  which  is  imposed  by  practi¬ 
cal  limitation  of  the  IP  solver,  the  ILO  based  algorithm 
produces  globally  optimal  solutions  with  very  high  fre¬ 
quency  (>  70%),  and  average  performance  is  within 
1.12%  of  the  optimal  solution. 

1.  Introduction 

Since  wireless  ad  hoc  and  sensor  networks  oper¬ 
ate  over  a  battery  energy  limited  environment,  the  de¬ 
sign  of  energy-efficient  network  algorithms  has  been  an 
active  research  topic.  Especially,  the  Minimum  En¬ 
ergy  Broadcast,  first  proposed  and  investigated  by 
Wieselthier  et  al.  [26,  27],  and  its  variant  Minimum 
Range  Assignment  [4,  5]  have  recently  attracted 
significant  attention  from  research  community,  since 
broadcast  is  a  major  communication  mode  for  informa¬ 
tion  dissemination.  Both  the  general  graph  and  geomet¬ 
ric  graph  version  of  this  problem  have  been  proven  to  be 
NP-complete  [3, 18].  Therefore,  unless  P  =  NP,  it  is 
not  likely  there  exist  exact  polynomial  time  algorithms 
to  solve  this  problem. 
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As  an  exact  method  to  obtain  globally  optimal  so¬ 
lutions,  a  mixed  integer  programming  (IP)  formulation 
for  Minimum  Energy  Broadcast  has  been  pre¬ 
sented  by  Das,  et  al.  [7].  Albeit  valuable  for  theo¬ 
retical  reasons,  its  practical  usage  is  limited  to  only 
small  input  sizes,  as  is  the  usual  case  for  many  IP  prob¬ 
lems.  To  deal  with  practical  input  size  we  still  need 
to  rely  on  efficient  heuristic  approximate  algorithms, 
which  usually  produce  satisfactory  results  in  polyno¬ 
mial  time  bound  but  without  the  guarantee  of  global 
optimality.  Over  the  past  years,  several  (greedy)  con¬ 
struction  heuristics  have  been  developed  such  that  a 
feasible  solution  is  built  from  scratch,  say,  using  lo¬ 
cations  of  nodes  as  input.  Some  of  the  representative 
algorithms  in  this  category  are:  the  Broadcast  Incre¬ 
mental  Power  (BIP)  [26],  Embedded  Wireless  Multi¬ 
cast  Advantage  (EWMA)  [3],  Broadcast  Average  In¬ 
cremental  Power  (BAIP)  [24],  Greedy  Perimeter  Broad¬ 
cast  Efficiency  (GPBE)  [13],  Center  Oriented  Broadcast 
Routing  Algorithm  (COBRA)  [16],  Steiner  tree  based 
heuristics  by  Liang  [18],  Ant  Colony  System  (ACS) 
based  approach  [6],  and  G-Remit  [25],  etc. 

In  the  general  field  of  combinatorial  optimization, 
the  iterative  improvement  heuristics  have  been  another 
class  of  approximate  algorithms,  where  an  initial  feasi¬ 
ble  solution,  possibly  obtained  from  a  construction  al¬ 
gorithm  or  generated  randomly,  is  iteratively  improved. 
Local  search  and  other  meta-heuristics  (e.g.,  simulated 
annealing  and  tabu  search)  belong  to  the  class  of  iter¬ 
ative  improvement  heuristics.  Many  variants  of  local 
search  are  known  to  be  state-of-the-art  heuristics  for  a 
variety  of  combinatorial  optimization  problems;  for  in¬ 
stance,  for  Traveling  Salesman  Problem  (TSP), 
while  there  exist  many  tour  construction  algorithms 
(e.g.,  Christofides’  heuristics,  see  [11]  for  a  survey),  the 
well-known  basic  fc-exchange  and  Lin-Kernighan  (LK) 
local  search  heuristics  [19,20],  developed  over  three 
decades  ago,  still  remain  to  be  the  algorithms  of  choice 
in  terms  of  efficiency  and  effectiveness. 

On  the  other  hand,  we  observe  however,  much  less 
attention  has  been  given  in  this  direction  for  Minimum 
Energy  Broadcast.  In  fact,  there  are  a  few  known 
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algorithms  that  may  be  considered  as  local  search  al¬ 
gorithms  in  a  broader  sense,  including  Post  Sweep  pro¬ 
cedure  [26],  Iterative  Maximum  Branch  Minimization 
(IMBM)  [17],  r-Shrink  [6],  although  not  specifically 
stated  as  such.  Interpreting  these  algorithms  as  local 
search  facilitates  analysis  of  advantages  or  disadvan¬ 
tages  therein,  and  after  the  underlying  concepts  of  each 
algorithm  are  distilled  as  local  search,  we  can  easily  ex¬ 
tend  the  concepts  to  devise  better  performing  heuris¬ 
tics.  Performance-wise,  above  heuristics  (IMBM,  r- 
Shrink,  Post-Sweep)  are  not  impressive,  resulting  in  lo¬ 
cal  optima  only  comparable  to  tree  construction  heuris¬ 
tics  (e.g.,  EWMA)  in  their  solution  quality. 

In  our  prior  work  [15],  we  presented  a  highly  effec¬ 
tive  heuristic  algorithm  (comparable  to  the  LK  search 
for  TSP  [20])  called  the  Largest  Expanding  Sweep 
Search  (LESS),  based  on  local  search  principle.  In  this 
paper,  we  extend  this  concept  and  present  an  iterated  lo¬ 
cal  optimization  based  heuristics  relying  on  two  neigh¬ 
borhood  structures,  which  are  designed  to  be  more  gen¬ 
eral  and  to  maximize  the  correlation  with  the  objective 
function  (the  meaning  of  which  will  be  clarified  later), 
still  maintaining  computational  efficiency.  One  short¬ 
coming  of  these  approaches  is  that  they  are  centralized. 
Unless  the  broadcast  routing  tree  is  computed  at  one  lo¬ 
cation  and  distributed  to  all  other  nodes,  their  practical 
usage  as  a  broadcast  routing  protocol  may  be  limited. 
However,  we  believe  the  quest  to  develop  ever  better 
performing  algorithms  has  never  ended  (as  is  the  case 
for  TSP  problem).  Eventually,  the  insights  gained  dur¬ 
ing  the  course  of  our  investigation  as  to  which  mecha¬ 
nism  provides  trees  yielding  smaller  cost  may  be  bene¬ 
ficial  in  other  related  research  activities,  including  dis¬ 
tributed  and  localized  implementations. 

The  remainder  of  this  paper  is  organized  as  follows. 
In  the  next  section,  a  system  model  used  throughout  this 
paper  and  the  problem  formulation  will  be  presented.  In 
Section  3,  we  briefly  review  the  concept  of  iterated  local 
optimization.  In  Section  4,  we  discuss  two  neighbor¬ 
hood  structures  suitable  for  MEB  problem,  which  are 
one  of  the  most  crucial  ingredients  in  iterated  local  op¬ 
timization.  Subsequently,  the  details  of  search  strate¬ 
gies  and  the  choice  of  initial  feasible  solution  are  dis¬ 
cussed  in  Section  5.  Simulation  results  and  conclusions 
are  provided  in  Section  6  and  7. 

2.  System  Model  and  Problem  Formulation 


quired  power  to  support  a  link  {u,  v)  from  node  u  to  v 
is  denoted  by  P^v ,  which  consists  of  RE  transmit  power 
and  signal  processing  power  for  signal  transmission  and 
reception.  Due  to  path  loss  in  a  wireless  channel,  the  re¬ 
ceived  power  level  is  inversely  proportional  to  distance 
and  attenuates  as  ||t6v||~“,  where  a  denotes  the  path 
loss  factor  which  usually  satisfies  2  <  a  <  7.  There¬ 
fore,  Puv  =  7  II Ml’ IP  +  A  is  a  relatively  accurate  link 
model  assuming  stationary  nodes  and  a  time-invariant 
(or  slowly  varying)  channel,  where  7  denotes  the  re¬ 
ceiver  sensitivity  threshold  and  A  is  a  constant  power 
consumption  factor  due  to  signal  processing  at  a  node. 
To  avoid  undue  complication  of  notations,  we  tenta¬ 
tively  assume  7=1  and  A  =  0.  We  note  that  our  al¬ 
gorithm  works  independent  of  an  individual  link  model. 
The  power  consumption  at  node  u  is  denoted  by  p  (u), 
and  its  corresponding  transmission  range  is  r  (m)  ,  which 
are  bounded  by  p  (u)  <  p^ax  and  r  (u)  <  r max,  respec¬ 
tively.  Therefore,  the  required  transmit  power  of  node 
u  to  reach  node  is  p  (u)  =  Puv 

Given  a  network  represented  by  a  directed  weighted 
graph  G  =  {V,  A,  c)  with  a  set  U  =  {1, 2, . . . ,  n}  of 
nodes  with  cardinality  n  =\V\  and  a  set  A  of  m  =  |7l| 
directed  edges  (links),  and  a  real  valued  cost  function  c  : 
A  ^  defined  on  every  edges  e  &  A,  the  Minimum 
Energy  Broadcast  (MEB)  problem  seeks  to  find 
a  s-arborescence  (served  as  a  broadcast  routing  tree) 
T  =  (y,  At)  rooted  at  source  node  s  that  can  reach 
every  node  v  G  U  in  a  wireless  network  either  directly 
or  through  intermediate  nodes  (in  multiple  hops)  with 
minimum  overall  power  cost.  Then  the  general  graph 
version  can  be  formulated  as: 


minimize 


E 

uGV 


max 

{u,v)gAt 


\/T  eTCG, 


(1) 


where  T  denotes  the  family  of  s-arborescences  (di¬ 
rected  trees)  rooted  at  s  and  spanning  all  nodes  in  V . 
The  geometric  version  of  MEB  assumes  the  graph  be  a 
geometric  graph;  That  is,  a  directed  edge  {u,v)  G 
exists  if  and  only  if  ||uv||  <  r  (u).  Note  (u,  v)  €  A  does 
not  necessarily  mean  (v,  u)  G  A,  since  r  (u)  ^  r  (v)  in 
general. 

In  what  follows,  we  will  use  the  following  notations: 


•  p  (u)  =  total  required  power  to  enable  communi¬ 
cations  between  nodes  lying  at  range  r  (u)  apart 
from  u. 


We  assume  that  each  node  (host)  in  a  wireless  ad  hoc 
network  is  equipped  with  an  omnidirectional  antenna, 
and  acquires  its  location  information  either  through 
GPS  or  other  localization  techniques  [2].  Let  ||mw||  be 
the  Euclidean  distance  between  node  u  and  v.  The  re¬ 


•  Puv  =  total  required  power  to  maintain  a  commu¬ 
nication  link  between  node  u  and  v. 

•  6t  (u)  =  {v  I  {u,v)  G  At},  the  logical  neighbor 
of  node  m  in  a  given  tree  T  which  the  set  of  ad¬ 
jacent  (child)  nodes  of  node  u.  The  cardinality  of 


(5t  (u)  corresponds  to  the  outdegree  of  node  u  in  S2  =LocalSearch(jF,  c,  A^2),  and  Ni  D  N2,  it  is 

the  tree  T,  i.e.,  deg  (tt)  =  I^t  (m)|  •  The  indegree  highly  likely  that  c(st)  <  c(s^),  although  it  is  not  al- 

of  every  node  is  always  1,  except  the  source  node  ways  guaranteed. 
s  for  which  the  indegree  is  0. 

•  ttt  (u)  =  {t;  I  (t;,M)  G  At}  ,Vu  G  y\{s}  and 
ttt  (s)  =  <t>,  the  unique  parent  node  of  node 
u  in  T.  We  also  use  the  notation  ttt  (S')  = 

•  A/’m  (v)  =  {w  I  0  <  ||mw||  <  lluvll  ,  w  G  yj ,  the 
physical  neighbor  of  node  u  that  is  the  set  of  all 
nodes  within  the  range  r  (u)  =  ||uv|| . 

•  =  the  set  of  nodes  lying  on  a  directed  path 
from  the  source  s  to  node  u,  where  s,  u  G 

3.  Iterated  Local  Optimization  (ILO) 

In  this  section,  a  brief  background  on  iterated  lo-  algorithms  rely  on  a  single  neigh- 

cal  optimization  is  presented.  Excellent  references  on  structure.  As  a  way  to  improve  upon  local 

basic  local  search  can  be  found  m  [1 1, 22].  For  more  to  combine  multiple 

recent  developments  m  meta-heunstics,  readers  are  re-  neighborhoods  for  descent;  e.g.  Martin  et  al.  [21]  pro- 

ferred  to  [8].  Given  an  optimization  problem  with  m-  ^  simulated  annealing  based  algorithm  for  TSP 

stances  {T,  c)  where  T  is  the  domain  of  feasible  points;  double-bridge  4-exhange  neighborhood  combined 

c:  :F^R  is  the  cost  function,  a  neighborhood  is  a  standard  local  search  using  3-opt  neighborhood, 

mapping  N  :  T 2^  defined  for  every  instance  in  ^  commonly  called  under  the  various 

Then,  a  feasible  solution  x  G  is  called  locally  [1 large-step  Markov 

optimal  with  respect  to  N  (or  simply  locally  optimal  pi],  chained  local  search  [1],  multilevel  local 

whenever  N  is  understood  by  context),  if  c{x)  <  c  (y) ,  [23]^  ^^^i^ble  neighborhood  search  [8].  In 

WyGN  (x),  and  is  called  globally  optimal,  if  c  (x)  <  p^p^^^  approach  under  the  more 

c  (y) ,  Vy  G  A-.  The  ‘standard’  local  search  algorithm  optimization.  For  in- 

for  a  minimization  problem  can  be  written  as:  -pSP  [1],  Iterated  FK 

(IFK)  [11],  Fin-Kernighan-Helsgaun  (LKH)  [9],  and 
multilevel  LKH  [23]  algorithms  are  based  on  the  prin¬ 
ciple  of  iterated  local  optimization,  and  constitute  the 
current  state-of-the-art  algorithms  in  terms  of  efficiency 
and  effectiveness. 

A  high  level  description  of  iterated  local  opti¬ 
mization  is  presented  in  Fig.  2,  which  utilizes  two 
neighborhood  structures  and  A^2-  First,  starting 
from  an  initial  feasible  solution  sq,  after  applying  Lo- 
calSearch(jF,  c,  A^i)  using  the  neighborhood  Ni,  we 
arrive  at  a  local  optimum  s*.  Next  we  apply  perturba- 
which  shows  that,  after  an  initial  feasible  solution  is  tion  or  ‘kick’  to  s*  using  Perturbation(s*,  A^2)  by  the 

generated  by  a  predefined  feasible  solution  genera-  neighborhood  N2  to  escape  the  local  optimum  (by  al- 

tor  GeneratelnitialFeasibleSolution,  this  solution  is  lowing  temporary  hill  climbing),  whose  outcome  is  an 

iteratively  refined  until  there  is  no  further  improve-  intermediate  stage  s'.  LocalSearch(jF,  c,  A^i)  is  once 

ment  within  the  neighborhood.  The  output  of  this  again  applied  but  now  using  s'  as  an  initial  feasible  so- 

procedure  is  a  local  optimum  s*,  which  is  heav-  lution,  and  accept  or  reject  the  output  depending  on  the 

ily  dependent  on  the  neighborhood  N  and  the  ini-  AcceptanceCriterion.  This  step  is  repeated  until  a  ter- 
tial  feasible  solution.  In  general,  the  larger  the  mination  condition  is  satisfied.  We  will  discuss  individ- 

neighborhood  size,  the  better  the  quality  of  the  lo-  ual  component  in  the  context  of  Minimum  Energy 

cal  optima,  i.e.,  if  s}  =LocalSearch(jF,  c,  iVi),  Broadcast  in  Section  5. 


Figure  1.  Pseudo-code  of  a  standard  local 
search  procedure. 


procedure /teratec? Local  Optimization{!F ,  c,  Ni,  N2) 

So  < —  GeneratelnitialFeasibleSolution 
s*  < —  LocalSearch(J^,  c,  iVi) 

repeat 

s'  < —  Perturbation(s*,  iV2) 
s*'  < —  LocalSearch(s',  c,  A^i) 
s*  < —  AcceptanceCriterion(s*,  s*') 
until  TerminationCondition  satisfied 
return  s* 

Figure  2.  Pseudo-code  of  an  iterated  local 
optimization  procedure. 


4.  Two  Neighborhood  Structures  for  Iter¬ 
ated  Local  Optimization 

The  local  optimization  for  Minimum  Energy 
Broadcast  (MEB)  requires  an  arborescence  as  an 
initial  feasible  solution  and  the  intermediate  improved 
solution  at  each  iteration.  To  preserve  tree  structure, 
the  basic  premise  is  that  every  node  should  be  con¬ 
nected  from  the  source,  and  there  should  be  no  cy¬ 
cle.  Given  a  tree  T,  we  can  uniquely  and  easily  iden¬ 
tify  the  corresponding  range  assignment  vector  r  = 
(r  (1) ,  r  (2) , . . . ,  r  (n))  G  TZ  C  M"  in  0  (n)  time, 
where  TZ  is  the  domain  of  feasible  solutions  for  Mini¬ 
mum  Range  Assignment  (MRA).  But  the  converse 
is  not  true,  because  a  group  of  trees  are  mapped  into  the 
same  range  assignment  vector  r.  See  Eig.  3  for  example. 
In  effect,  the  domain  of  feasible  solutions  TZ  for  MRA 
is  a  partition  of  the  domain  of  the  feasible  solutions  T 
for  MEB.  Thus  if  we  deal  with  range  assignments  f  in¬ 
stead  of  trees  T,  the  search  space  of  local  search  [1 1,22] 
for  MEB  problem  may  be  significantly  simplified. 


Figure  3.  Range  assignment  vs.  tree. 

(a)  A  tree  and  its  range  assignment,  (b) 
For  the  same  range  assignment,  muiti- 
pie  trees  are  aiiowabie.  For  this  specific 
topoiogy,  2*  =  256  trees  are  mapped  to 
the  same  range  assignment  r{s)  —  ||su||, 
r{u)  —  ||itu||.  Unaffected  edges  are  not 
dispiayed. 

In  general,  it  is  suggested  that  neighborhoods  be  de¬ 
signed  in  order  to  induce  maximum  correlation  between 
objective  function  values  of  adjacent  points  [11].  In  or¬ 
der  to  maximize  correlation  between  objective  function 
and  neighborhood  function,  it  is  better  to  directly  deal 
with  range  assignments  instead  of  trees.  As  observed  by 
Wieselthier  et  al.  [26],  this  also  matches  well  with  the 
general  philosophy  that  node-based  approach  (range)  is 
preferred  to  link-based  approach  (tree)  in  wireless  envi¬ 
ronment.  This  observation  is  the  basis  of  our  new  ap¬ 
proach  and  leads  to  an  algorithm  that  converges  quickly 
to  a  high  quality  local  optimum.  In  the  following,  if  at 
most  k  ranges  of  nodes  are  allowed  to  increase  from  a 


given  range  assignment  r,  we  call  it  an  order-fc  (range- 
based)  neighborhood  and  denote  hy  NH  (r),  i.e., 

(r) 

=  {r'  \  r'  (v)  >  r  (v)  ,yv  G  S  C  V  A  [S']  <  k}  . 

4.1.  Expanding  Sweep  Search  Neighborhood 

First  we  present  an  order- 1  range-based  neighbor¬ 
hood  structure  called  the  Expanding  Sweep  Search 
{ESS)  neighborhood  which  also  preserves  tree  connec¬ 
tivity  from  source  s.  The  idea  of  ESS  is  simple;  in 
exchange  for  an  increase  in  a  node’s  power  (or  equiv¬ 
alently  range),  how  much  reduction  in  other  nodes’ 
power  can  be  achieved  is  tested.  Before  we  proceed 
further,  we  note  that  all  operations  in  our  algorithm  are 
based  on  the  following  trivial  inequality:  for  any  finite 
sets  of  real  numbers  and  S2  such  that  C  S2, 

maxis'll  <  max {S'2}  .  (2) 

Thus  it  works  on  general  graphs  with  arbitrary  edge 
costs,  even  in  the  case  of  asymmetric  edge  costs.  How¬ 
ever,  if  a  geometric  graph  model  is  used,  we  can  lever¬ 
age  its  geometric  properties  to  find  more  efficient  algo¬ 
rithms. 

The  ESS  transformation  may  be  best  explained  with 
the  illustrations  presented  in  Fig.  4.  Fig.  4(a)  shows  the 
original  tree  T,  where  paths  from  source  s  to  nodes  c, 
e,  /  are  omitted.  The  direction  of  each  edge  is  always 
assumed  away  from  s.  The  transmission  range  of  each 
node  is  also  drawn.  Let  us  keep  focus  on  node  u,  which 
originally  has  the  range  r  (m)  =  ||u(i||.  Suppose  now 
node  u  decides  to  expand  its  range  to  reach  a  destination 
node  V  such  that  ||mi;||  >  ||Mfi||,  at  the  increase  in  total 
cost  by 

A+  (u)  =  -  Pw  >  0  (3) 

which  is  called  the  incremental  power.  The  expanded 
range  is  illustrated  as  a  solid  line  circle  in  Fig.  4(b). 
Then  every  node  within  the  range,  i.e.,  the  physical 
neighbor  A/'u  {v)  =  {a,b,d,g,h,  i,j,  u,  r;},  can  hear  the 
broadcast  message  from  s  through  node  u.  Thus,  re¬ 
gardless  of  previous  parent  node,  each  node  in  Afu  {v) 
may  update  its  parent  node  to  be  node  u,  since  they  can 
get  the  message  from  node  u.  However,  cautions  should 
be  given  since  if  any  of  the  path  nodes  from  the  source 
to  node  u,  denoted  by  Hg^^i  =  {s, . . .  ,c,b,a,u}, 
would  update  its  parent  to  node  u,  the  path  from  s  to 
u  would  be  broken  and  the  tree  might  become  discon¬ 
nected.  Therefore,  only  the  nodes  {d,  g,  h,  i,j,  u},  the 
physical  neighbors  other  than  the  nodes  in  Hg^^i,  are 
allowed  to  become  the  children  of  node  u.  That  is,  ev¬ 
ery  node  in  A4„  (v)  =  Afu  (v)  iHg^^  (marked  with  an 


Figure  4.  Illustration  for  neighborhood  structures  (a)  Original  tree  T.  (b)  Node  u  expands 
its  range  to  reach  node  v.  (c)  A  neighbor  tree  T'  of  T  in  Expanding  Sweep  Search  (ESS) 
neighborhood. 


empty  circle  in  Fig.  4)  updates  its  parent  to  node  u.  The 
resulting  transformed  tree  T'  by  this  parent  updates  is 
shown  in  Fig.  4(c). 

Now  let  us  analyze  how  much  reduction  in  total 
tree  cost  we  can  achieve  by  this  transformation.  Note 
that  updating  parents  of  node  in  (n)  in  turn  gives 
a  chance  for  the  original  parents  w  G  Qu  (f)  — 
7r(Alii  (f  ))\  {rt}  =  {h,  d,  e,  /}  marked  with  square  in 
Fig.  4  to  reduce  their  ranges  resulting  in  power  reduc¬ 
tion  by 

A~{w)=p{w)-  max  {Pwk}  (4) 

fcG(5T  u  (l^) 

=  max  {Pnik}  -  max  {Pwk}  >  0, 

k^ST{w)  k^ST{w)\A4u{'^) 

since  6t  (w)  \Mu  {v)  C  6t  (w)  and  using  (2).  Ac¬ 
counting  for  every  affected  node  w  G  Qu{v),  the  pure 
gain  Q  {u,  v)  in  total  cost  is 

g  {u,  v)  =  ^  A"  (w)  -  A+  (u) ,  (5) 

weQu{v) 

where  the  first  term  on  the  right  hand  side  is  called 
sweeping  gain.  We  denote  the  ESS  transformation  by 
the  postfix  form  T'  =  T  \  \  £SS  {u,v).  Then  the  ex¬ 
panding  sweep  search  neighborhood  denoted  by  N Hess 
is  defined  as 

NHess  (T)  =  {T'\r  =  T\\  £SS  {u,  v)  ,Wu,veV}. 

Since  \NHess  (T)|  <  (2),  the  size  of  neighborhood 
is  O  (n^)-  After  the  expanding  sweep  search  neigh¬ 
borhood  NHess  is  clearly  defined,  the  LESS  algorithm 
[15]  is  nothing  more  than  LocalSearch(T,  c,  NHess) 
with  the  steepest  descent  search  strategy.  That  is,  given 
an  initial  feasible  solution,  the  best  neighbor  among 
NHess  is  chosen  at  each  iteration,  until  there  is  no  more 
gain  in  cost. 


procedure  Expanding  Sweep  Search  (ESS)  (u,  v) 
Input:  a  s-arborescence  T  G  T 

1.  T'  ^  T 

2.  expand  range  of  node  u  to  n,  i.e.,  r(u)  <—  ||Mn||, 
Vv  G  V  s.t.  ||■u^;||  >  r{u),v  ^  Ifs^^ 

3.  Mu  {v)  ^  Afu  (v)  \ns^u 

4.  TTt'  (Mu  (v))  ^  u 

5.  return  T' 


Figure  5.  Pseudo-code  for  Expanding 
Sweep  Search  (ESS)  transformation. 


4.2.  Edge  Exchange  Neighborhood 

An  alternative  (and  traditional)  approach  is  to  con¬ 
sider  the  tree  structure  itself  instead  of  range  assign¬ 
ments.  The  most  natural  neighborhood  structure  con¬ 
ceivable  when  the  tree  structure  is  a  feasible  solution  is 
the  elementary  tree  transformation  for  undirected  graph 
as  described  in  Papadimitriou  and  Steiglitz  (see  Chap¬ 
ter  19  in  [22]);  that  is,  add  a  non-tree  edge  to  a  tree 
and  remove  an  edge  from  the  unique  cycle  formed  by 
previous  edge  addition.  Because  we  are  dealing  with  a 
rooted  directed  tree  (arborescence)  for  Minimum  En¬ 
ergy  Broadcast,  we  can  only  apply  elementary  tree 
transformation  on  the  underlying  graph  of  the  directed 
graph.  Therefore  we  call  it  an  Edge  Exchange  Neigh¬ 
borhood  (NHee)  that  is  an  adaptation  of  the  elementary 
tree  transformation  to  directed  graphs. 

An  equivalent  interpretation  of  edge  exchange  neigh¬ 
borhood  is  as  follows:  Remove  an  edge  (u,  v)  from 
T  which  partitions  the  tree  into  two  subtrees  Ti  = 
(Suv,  Atx  )  and  T2  =  (V\Suv,  )  where  Suv  denotes 
the  connected  component  containing  s  and  u.  Adding 
an  edge  (u',  v')  G  Suv  x  V\Suv  to  Ti  ur2  and  updating 


the  direction  of  edges  in  T2  forms  a  neighbor  of  T: 

NHee  (T)  =  {T'  I  T'  =  {V,  At\  {u,  v)  U  {u',  v'))  , 
V  (u,v)  G  At,  V  {u',v')  G  Suv  x  V\Suv, 

reverse  edge  direction  }. 

The  size  of  neighborhood  is  (|V^|  —  |5'„t,|)  = 

O(n^).  Checking  the  cost  of  a  tree  c(T),  using  a 
simple  array  as  the  data  structure  for  a  tree,  requires 
O  (n)  running  time.  Hence  if  we  use  the  algorithm 
LocalSearch(T,  c,  A^iTee)  with  steepest  descent 
strategy,  each  iteration  may  take  O  (n^),  which  is 
quite  expensive.  In  the  course  of  updating  the  direction 
of  edges  in  T2,  it  is  possible  that  O  (n)  nodes  may 
be  required  to  increase  their  ranges.  Therefore,  this 
neighborhood  is  complementary  to  NHess-  Other  pos¬ 
sible  extension  to  NHee  neighborhood  is  to  combine 
NHee  and  NHess-  that  is,  for  every  expanded  range 
in  NHee,  we  can  apply  the  ESS  transformation. 
Although  the  combined  neighborhood  may  give  better 
performance,  running  time  per  iteration  is  costly. 
Trade-offs  in  neighborhood  size,  time  complexity  and 
the  quality  of  local  optima  is  required.  Therefore, 
we  restrict  the  usage  of  NHee  as  a  perturbation 
method  for  iterated  local  optimization.  Since  Lo- 
calSearch(T,  c,  NHess)  7^LocalSearch(T,  c,  NHee) 
in  general,  they  are  hence  mutually  beneficial  in  en¬ 
hancing  the  quality  of  local  optimum  obtained  by  the 
other. 

Theorem  1  The  following  properties  of  these  two 
neighborhoods  satisfy: 

(a)  Expanding  Sweep  Search  neighborhood  encom¬ 
passes  neighborhood  for  Post  Sweep  algorithm. 

(b)  Edge  Exchange  neighborhood  encompasses  neigh¬ 
borhoods  for  IMBM  and  1 -Shrink. 

Proof,  (a)  Without  range  expansion  in  NHess,  inspect¬ 
ing  gain  for  current  range  assignment  is  no  different 
from  the  post  sweep  algorithm  [26].  (b)  The  neighbor¬ 
hoods  used  in  Iterative  Maximum  Branch  Minimization 
(IMBM)  [17]  and  r-shrink  algorithms  [6]  are  obtained 
by  updating  a  node’s  parent  to  a  node  which  is  not  one 
of  its  descendants.  Clearly,  this  is  a  special  case  of 
NHee  when  v  =  v'.  ■ 

Therefore,  every  local  optimum  by  NHess  (or  NHee) 
is  also  a  local  optimum  by  Post  Sweep  (or  IMBM  and  1- 
Shrink),  respectively.  Furthermore,  note  that  although 
T'  in  Fig.  4(c)  is  a  one-hop  neighbor  of  T  in  NHess, 
the  corresponding  tree  structure  is  large  number  of  hops 
away  in  terms  of  edge  exchange  neighborhood  NHee- 
This  is  the  true  advantage  of  using  range  assignments 


in  our  algorithms,  because  the  number  of  required  iter¬ 
ation  can  be  significantly  reduced,  compared  to  using 
tree  structure,  without  much  added  complexity. 

5.  Search  Strategies  and  Algorithm  Details 

Our  current  implementation  of  iterated  local  opti¬ 
mization  in  this  work  strictly  follows  the  reference  high 
level  description  in  Fig.  2.  It  is  based  on  deterministic 
local  search  LocalSearch(T,  c,  NHess),  and  a  random 
perturbation  Perturbation(s*,  NHee)-  We  will  discuss 
each  component  in  Fig.  2  in  the  order  of  appearance. 
Before  we  further  proceed,  let  us  discuss  some  of  the 
underlying  philosophies  of  our  algorithm. 

The  underlying  philosophy  of  our  current  work  is 
exactly  the  opposite  of  the  Iterative  Maximum  Branch 
Minimization  (IMBM)  [17]  and  r-shrink  [6]  algo¬ 
rithms.  As  suggested  by  the  names  of  these  algorithms, 
“range  reduction”  is  emphasized,  meaning  that  the  max¬ 
imum  branch  (edge)  and  the  corresponding  range  of 
each  node  is  progressively  broken  into  smaller  pieces 
(i.e.,  into  multihop  edges  or  ranges).  For  instance  in 
IMBM,  as  an  initial  feasible  solution,  a  single  maxi¬ 
mum  range  from  the  source  s  covering  the  rest  of  the 
nodes,  i.e.,  r  (s)  =  max„gy  {||sr||}  and  the  corre¬ 
sponding  star  topology  rooted  at  s,  is  used.  It  is  it¬ 
eratively  reduced  and  transformed  to  different  multi¬ 
hop  tree  structure  [17].  The  r-shrink  algorithm  resem¬ 
bles  the  IMBM,  but  relies  on  a  deterministic  maximum 
branch  reduction  method. 

On  the  other  hand,  our  LESS  algorithm  [15]  progres¬ 
sively  expands  the  ranges  of  nodes.  While  range  ex¬ 
pansion  and  range  reduction  are  like  the  flip  side  of  the 
coin  (one  can  be  achieved  at  the  cost  of  the  other),  the 
implication  in  local  search  is  radically  different.  That 
is,  to  reduce  the  range  of  a  node,  we  need  to  specifi¬ 
cally  search  for  the  nodes  whose  increase  in  range  can 
reduce  the  total  cost.  However,  as  observed  in  expand¬ 
ing  sweep  search,  expansion  in  a  node’s  range  naturally 
determines  which  nodes  can  reduce  their  ranges  and 
do  not  require  further  computations,  which  is  one  ad¬ 
vantage  of  our  approach.  Furthermore,  there  are  other 
reasons  range  expansion  may  be  generally  preferable  to 
range  reduction  for  MEB. 

Let  us  consider  the  topology  illustrated  in  Fig.  6, 
where  two  groups  of  m  nodes  lie  at  the  inner  and  outer 
circle  of  radius  ri  and  r2,  respectively,  where  |r2  —  ri  | 
is  assumed  very  small.  Suppose  the  configuration  in 
Fig.  6(a)  is  given  as  the  initial  feasible  solution,  as  is  the 
case  for  IMBM.  While  the  net  difference  between  two 
tree  cost  may  be  small  (r^  vs.  r“  -|-  m  (r2  —  ri)“),  ei¬ 
ther  deterministically  or  probabilistically  searching  for 
the  occasions  which  enable  reduction  in  range  r  (s)  is 


Figure  6.  Illustration  of  a  stable  constella¬ 
tion  demonstrating  that  reduction  in  range 
is  difficult  in  general. 


extremely  difficult.  That  is,  only  if  all  m  interior  nodes 
simultaneously  increase  their  ranges,  we  can  jump  out 
of  the  local  optimum  to  convert  to  the  better  configu¬ 
ration  in  Fig.  6(b).  In  terms  of  the  number  of  hops 
in  neighborhood,  these  two  solutions  are  far  apart  and 
hence  we  may  assume  there  exists  a  high  potential  bar¬ 
rier  in  between  these  configurations.  In  the  course  of 
range  reduction,  similar  situation  occurs  very  frequently 
and  hence  it  is  easy  to  fall  into  a  poor  local  optimum. 

GeneratelnitialFeasibleSolution:  Now  it  becomes 
quite  clear  which  initial  feasible  solution  is  a  suitable 
choice.  Note  that  the  stable  configuration  discussed 
above  may  have  been  avoided  from  the  outset  (although 
not  completely),  if  we  have  chosen  MST  as  the  initial 
solution.  In  addition,  since  MST  consists  of  short  length 
edges  [12],  it  provides  richer  search  space  for  range  ex¬ 
pansion  algorithms  as  ours.  Also  approximation  ratio  of 
MST  is  known  to  be  12  [24].  By  improving  upon  MST, 
the  same  worst  case  bound  is  guaranteed.  For  the  same 
reason,  both  BIP  and  SPT  are  equally  viable  choices 
but  not  considered  in  this  paper.  Certainly,  repeatedly 
restarting  from  a  randomly  generated  tree  as  a  starting 
point  can  only  make  our  results  stronger. 

LocalSearch:  As  the  LocalSearch  algorithm  in 

Fig.  2,  we  could  use  our  LESS  algorithm  [15].  How¬ 
ever,  we  can  speed  up  convergence  significantly  by  se¬ 
lecting  the  first  k  largest  gain  at  each  iteration  as  long  as 
the  ranges  are  disjoint,  for  which  we  call  the  variable- 
order  disjoint  (VOD)  steepest  descent  search  strategy. 
This  strategy  in  some  sense  resembles  the  variable- 
depth  search  by  Lin  and  Kernighan  [20]  in  that  indef¬ 
inite  orders  are  considered  as  long  as  the  gain  is  kept 
positive.  We  will  call  this  algorithm  as  LESS„od-  The 
details  are  listed  in  Fig.  7. 

Perturbation:  As  discussed  in  Section  4.2,  we  adopt 
edge  exchange  neighborhood  A^iFee  as  our  perturbation 


procedure  Variable-Order  Disjoint  LESS  (LESSvod) 

Input:  an  s-arborescence  t 

/*  requires  searching  whole  neighborhood  */ 

1.  fi{u) -f—  argmax  {Q{u,v)},'iu  €  V 

{v:||ui!||>r(n),t?^ns^.u} 

2.  sort  ui,  U2,  ■  •  ■ ,  Uk,  k  <  n  so  that 

g  {ui,  ^  (ui))  >  g  {uj,  ^  {uj))  >  0, 

\/l  <  i  <  j  <  k 

5.  R  < —  {ui} 

6.  for  i  =  2, ...,  k 

/*  if  ranges  are  disjoint  */ 

7.  if  ||uiu||  >  \\uifi{u^)\\  ||u^(u)||,Vu  G  R 

8.  R  < —  R  U  {ui} 

9.  end 

/*  iteratively  apply  ESS  |i?|  times  */ 

10.  t'  < —  t  II  ESS  {v,  (v)),  Vu  G  i? 

11.  if  c  {t')  <  c{t) 

/*  if  there  is  gain,  repeat  iteration  */ 

12.  t  < —  t'\  goto  step  1 

13.  return  t' 


Figure  7.  Pseudo-code  for  LESSvod  algo¬ 
rithm. 


method.  Due  to  its  O  (n^)  size  neighborhood,  steep¬ 
est  descent  type  search  strategy  can  be  very  expensive. 
Therefore,  random  perturbation  in  NHee  is  used  in  this 
paper,  which  requires  only  0(1)  running  time.  This  al¬ 
lows  most  of  the  time  in  iterated  local  optimization  will 
be  spent  by  LocalSearch  algorithms. 

AcceptanceCriterion:  As  in  simulated  annealing  or 
large-step  Markov  chain  [21],  the  following  acceptance 
criterion  is  popularly  used: 

AcceptanceCriterion  (s*,  s*') 
f  s*'  if  c  (s*')  <  c(s*) 

I  c(s*)— c(s*^) 

=  ■1  s*'  if  c  (s*')  >  c  (s*) ,  w.p.  e  T' 

I  s*  otherwise 

where  deteriorating  perturbed  intermediate  solution  s  *' 

is  accepted  with  probability  (w.  p.)  exp{ - — 1}, 

and  T  is  a  tunable  parameter  called  temperature.  If 
T  —  0,  only  the  better  solution  is  accepted.  If  T  is  very 
high,  random  walk  type  behavior  can  be  accomplished. 
In  this  paper,  we  adopt  this  acceptance  criterion  with 
T  —  0.  Simulations  using  T  >  0  and  changing  temper¬ 
ature  over  time  (cooling  schedule)  will  be  a  direct  but 
interesting  extension  to  current  work. 

TerminationCondition:  We  can  terminate  the  it- 


eration:  (i)  if  s*  has  reached  the  global  optimum, 
(ii)  if  maximum  pre-specified  number  of  iteration  has 
reached,  or  (iii)  if  pre-specified  time  has  elapsed.  In 
case  the  globally  optimal  solution  is  known  (say,  using 
an  IP  solver),  we  use  (i);  otherwise,  we  use  (iii)  with 
maximum  time  limit  set  to  60  seconds  using  a  Pentium 
1 .6GHz  equipped  machine.  Now  that  all  components  in 
iterated  local  optimization  is  established,  we  proceed  to 
simulations  results. 

6.  Simulation  Results 

Over  1000x1000  m^  deploy  region,  network  con¬ 
figurations  are  generated  by  uniformly  distributing  n 
nodes,  one  of  the  node  is  selected  as  the  source  s,  and 
broadcast  routing  trees  rooted  at  s  are  constructed.  By 
varying  the  number  of  nodes  within  the  deploy  region, 
we  have  effectively  changed  the  node  density  per  unit 
area.  The  path  loss  factors  a  =  2  and  4  are  used  in  our 
simulations.  For  each  network  instance,  globally  opti¬ 
mal  solutions  are  obtained  following  the  IP  formulation 
in  [7].  As  an  IP  solver,  CPLEX  [10],  one  of  the  most  so¬ 
phisticated  commercial  IP  solver,  was  used.  From  our 
experience,  a  network  size  of  up  to  20  nodes  is  the  prac¬ 
tical  maximum  limit  solvable  in  a  reasonable  time  (<  3 
hours  per  instance).  Flowever,  CPLEX  being  a  general 
purpose  IP  solver,  it  does  not  rule  out  the  possibility 
that  one  may  develop  customized  solver  by  generating 
problem-specific  cuts  for  branch-and-cut  methods  [22] 
either  to  speed  up  the  calculation  or  to  allow  for  larger 
network  sizes. 

For  small  network  sizes  4  <  n  <  20,  we  are  able  to 
compare  our  results  with  optimal  solutions  by  CPLEX, 
and  they  are  summarized  in  Table  1  and  2.  The  expres¬ 
sions  to  obtain  the  values  in  each  column  is  summarized 
underneath  the  tables.  Each  entry  in  the  tables  is  an  av¬ 
erage  of  1000  randomly  generated  network  instances. 
We  note  that  obtaining  the  whole  globally  optimal  solu¬ 
tions  took  several  weeks  using  a  Pentium  2.8GHz  pro¬ 
cessor  equipped  machine,  which  demonstrates  imprac¬ 
ticableness  of  the  IP  approach  in  real  situations.  All 
other  algorithms  are  implemented  in  Matlab  environ¬ 
ment.  We  allowed  up  to  60  seconds  for  the  iterated  local 
optimization  (ILO)  using  N Hess  +  NHee  to  run  before 
interrupted.  As  the  performance  comparison  measure, 
we  presented  both  the  average  total  transmit  power  and 
the  average  percentage  excess  (P.E.)  over  global  op¬ 
tima  defined  as 

^  ■=!  -  1]  X  100%, 

for  the  algorithm  algo.  Let  popt  denote  the  proba¬ 
bility  (frequency)  of  local  optima  by  ILO  being  actu¬ 


ally  global  optima  (OPT).  The  values  in  second  paren¬ 
thesis  in  column  (c)  roughly  represents  (l—popt)  = 
Pr  [c(Tilo)  >  c(Topt)]  in  percentage. 

We  can  observe  that  ILO  consistently  produces  very 
high  quality  local  optima,  among  which  globally  op¬ 
timal  solutions  are  obtained  with  high  frequency.  As 
noted  earlier,  NHee  is  complementary  to  N Hess ,  and 
hence  using  NHee  converts  a  large  fraction  of  subopti- 
mal  cases  into  global  optima.  Eor  instance,  for  a  =  2, 
n  <  20,  we  get  P.E.  <  1.117%  and  (1  —  Popt)  < 
30.8%.  In  case  of  a  =  4,  we  can  get  global  op¬ 
tima  for  almost  all  instances,  with  popt  >  96.5%  and 
P.E.  <  0.039%  for  n  <  20.  Out  of  interest  we  also 
tested  local  search  using  NHee  alone,  and  its  results 
are  reported  in  column  (d).  We  can  see  that  the  result  is 
not  impressive  enough,  although  still  better  than  IMBM 
or  1 -Shrink.  However,  we  should  note  that  this  level  of 
performance  is  achieved  without  random  restarts  as  in 
IMBM.  The  reason  we  get  better  performance  for  a  =  4 
is  quite  clear.  As  observed  in  [3, 14, 26],  the  larger  the 
path  loss  factor  a,  the  larger  the  penalty  for  increasing 
the  transmit  power.  As  the  wireless  broadcast  advantage 
becomes  negligible,  more  link-based  characteristics  is 
emphasized.  In  such  case  even  local  search  using  NHee 
alone  performs  reasonably  well  with  P.E.  <  1.887%, 
since  NHee  is  a  link-based  neighborhood.  Readers  may 
also  find  interest  in  the  performance  of  BIP  or  MST  rel¬ 
ative  to  the  global  optimum. 

Eor  a  larger  network  size  n  >  20,  comparison 
with  optimal  solutions  by  the  IP  solver  was  not  possi¬ 
ble.  Hence  as  a  reference,  performance  comparison  of 
LESS,  MST,  BIP  [26],  and  EWMA  [3]  are  presented 
in  Pig.  8.  Extrapolating  the  performance  level  of  LESS 
demonstrated  in  Table  1  and  2  to  larger  network  sizes, 
Pig.  8  provides  a  valuable  measure  how  other  algo¬ 
rithms  perform  relative  to  the  optimal.  Average  total 
transmit  power  is  compared  in  Pig.  8(a)  and  (c),  and 
the  ratio  of  them  is  compared  in  Pig.  8(b)  and  (d)  for 
a  —  2  and  4.  In  all  figures,  it  is  easy  to  notice  that 
LESS  consistently  outperforms  others  by  a  wide  mar¬ 
gin,  even  EWMA,  which  used  to  be  state-of-the-art. 
Eor  instance,  LESS  performs  approximately  25~31% 
(7.5~8.5%)  better  than  MST,  over  the  entire  range  from 
n  =  20  to  1000  for  a  =  2  (resp.  a  =  4).  Eor  a  =  4, 
due  to  high  dynamic  range,  it  is  hard  to  appreciate  the 
performance  gain  in  Fig.  8(c)  but  it  is  clearly  visible  in 
Fig.  8(d).  Clearly,  this  suggest  we  should  densely  de¬ 
ploy  network  nodes,  which  has  more  favorable  effect 
than  efficient  algorithms  for  large  a.  For  many  practi¬ 
cal  situations  when  a  high  quality  local  optimum  usually 
suffices,  our  work  can  be  used  for  benchmarking  other 
algorithms  for  comparison. 


Table  1.  Comparison  of  local  search  using  NHess  and  NHee  with  global  optima  {a  =  2) 

(a)n  (b)  OPT  (RE.)  (c)  ILO  (RE.)  (1  -  (d)  (P-E.)  (e)  BIP  (P.E.)  (Q  MST  (RE.) 


4 

4.3143e+5 

5 

4.3102e+5 

6 

4.2652e+5 

7 

4.2720e+5 

8 

4.2408e+5 

9 

4.1824e+5 

10 

4.1977e+5 

11 

4.1404e+5 

12 

4.1079e+5 

13 

4.0735e+5 

14 

4.0776e+5 

15 

4.0603e+5 

20 

3.9952e+5 

4.3143e+5  (0.0000%)  ( 0.0%) 
4.3102e+5  (0.0000%)  ( 0.0%) 
4.2655e+5  (0.0059%)  ( 0.2%) 
4.2734e+5  (0.0297%)  ( 0.7%) 
4.2444e+5  (0.0854%)  ( 2.0%) 
4.1850e+5  (0.0562%)  ( 2.0%) 
4.2045e+5  (0.1762%)  ( 3.9%) 
4.1513e+5  (0.2542%)  ( 6.1%) 
4.1195e+5  (0.2723%)  ( 6.8%) 
4.0853e+5  (0.2748%)  ( 7.3%) 
4.1012e+5  (0.5679%)  (13.8%) 
4.0853e+5  (0.6181%)  (15.5%) 
4.0390e+5  (1.1170%)  (30.8%) 


4.3455e+5  (0.7232%) 
4.3787e+5  (1.5893%) 
4.3691e+5  (2.4360%) 
4.4146e+5  (3.3380%) 
4.4421e+5  (4.7467%) 
4.4001e+5  (5.2051%) 
4.4510e+5  (6.0343%) 
4.4061e+5  (6.4173%) 
4.4078e+5  (7.3006%) 
4.4168e+5  (8.4276%) 
4.4337e+5  (8.7331%) 
4.4256e+5  (8.9969%) 
4.4340e+5  (10.9832%) 


4.4192e+5  ( 2.5994%) 
4.5459e+5  ( 5.9326%) 
4.5998e+5  (  8.1004%) 
4.7062e+5  (10.5279%) 
4.7761e+5  (13.3438%) 
4.7966e+5  (15.8574%) 
4.8863e+5  (17.2536%) 
4.8709e+5  (18.2900%) 
4.8602e+5  (19.1320%) 
4.8940e+5  (20.8748%) 
4.9478e+5  (22.1620%) 
4.8921e+5  (21.1684%) 
4.983 le+5  (25.4930%) 


4.7631e+5  (10.4685%) 
4.9352e+5  (15.5838%) 
5.0399e+5  (18.8144%) 
5.1693e+5  (22.3715%) 
5.2574e+5  (25.2631%) 
5.2605e+5  (27.6818%) 
5.3655e+5  (29.3122%) 
5.3244e+5  (29.6363%) 
5.335  le+5  (31.2965%) 
5.3479e+5  (32.5564%) 
5.3809e+5  (33.2930%) 
5.3572e+5  (33.1597%) 
5.3865e+5  (36.0863%) 


Table  2.  Comparison  of  local  search  using  NHess  and  NHee  with  global  optima  (a  =  4) 

(a)n  (b)  OPT  (P.E.)  (c)  ILO  (P.E.)  (1  -  p^pt)  jd)  NHee  (P-E.)  (e)  BIP  (P.E.)  (f)  MST  (P.E.) 


5 

1.3383e+ll 

6 

1.1127e+ll 

7 

1.0075e+ll 

8 

0.8822e+ll 

9 

0.781  le+11 

10 

0.7295e+ll 

11 

0.6441e+ll 

12 

0.5806e+ll 

13 

0.5431e+ll 

14 

0.4971e+ll 

15 

0.4571e+ll 

20 

0.3303e+ll 

1.3383e+ll  (0.0000%)  (0.0%) 
1.1127e+ll  (0.0000%)  (0.0%) 
1.0074e+ll  (0.0009%)  (0.1%) 
0.8825e+ll  (0.0397%)  (0.2%) 
0.7813e+ll  (0.0207%)  (0.4%) 
0.7297e+ll  (0.0412%)  (0.8%) 
0.6444e+!l  (0.0438%)  (0.9%) 
0.5808e+ll  (0.0384%)  (0.6%) 
0.5432e+ll  (0.0418%)  (1.5%) 
0.4972e+ll  (0.0468%)  (1.4%) 
0.4573e+ll  (0.0645%)  (2.1%) 
0.3305e+ll  (0.0391%)  (3.5%) 


1.3448e+ll  (0.4205%) 
1.1187e+ll  (0.5850%) 
1.0176e+ll  (0.8120%) 
0.8932e+ll  (1.1653%) 
0.7895e+ll  (0.9693%) 
0.7367e+ll  (1.0714%) 
0.6522e+ll  (1.2194%) 
0.5891e+ll  (1.3915%) 
0.5534e+ll  (1.7298%) 
0.5045e+ll  (1.4731%) 
0.4647e+ll  (1.7533%) 
0.3365e+ll  (1.8868%) 


1.3686e+ll  (2.8603%) 
1.1504e+ll  (3.5228%) 
1.0476e+ll  (4.2778%) 
0.9249e+ll  (4.9752%) 
0.8223e+ll  (5.8491%) 
0.7713e+ll  (6.0221%) 
0.6837e+ll  (6.2592%) 
0.6174e+ll  (6.3946%) 
0.5796e+ll  (6.8126%) 
0.5318e+ll  (7.2157%) 
0.4875e+ll  (6.8394%) 
0.3529e+ll  (7.2012%) 


1.4025e+ll  (5.3330%) 
1.1854e+ll  (6.4568%) 
1.0713e+ll  (7.0146%) 
0.9527e+ll  (8.3369%) 
0.8418e+ll  (8.7395%) 
0.7916e+ll  (9.1413%) 
0.7015e+ll  (9.0619%) 
0.6336e+ll  (9.3804%) 
0.5912e+ll  (9.1565%) 
0.5434e+ll  (9.9587%) 
0.4987e+ll  (9.5634%) 
0.3663e+ll  (11.2691%) 


(a)n 


network  size  (b) 


lE 


j  =  l  CjiToPj)  (C)  ^  Ej  =  l  Cj{TNHess+NHee)  (d)  \  E 
i  E)=i  Cj{Tbip)  (f)  J  E)=i  Cj{Tmst) 


k 

7  =  1 


Cj{TNHeJ  (e) 


7.  Conclusions 

We  investigated  efficient  heuristics  for  the  minimum 
energy  broadcast  (MEB)  problem.  While  there  have 
been  much  efforts  in  developing  efficient  tree  construc¬ 
tion  algorithms  beginning  from  the  BIP  algorithm,  we 
provided  different  algorithms  based  on  an  iterative  im¬ 
provement  approach.  Specifically,  fhrough  the  use  of  it¬ 
erated  local  optimization  techniques,  it  is  demonstrated 
that  the  performance  of  our  former  algorithm.  Largest 
Expanding  Sweep  Search  (LESS),  can  be  further  en¬ 
hanced.  Extensive  comparison  study  with  globally  op¬ 
timal  solutions  obtained  by  an  integer  programming 
solver  confirmed  fhat  the  produced  solutions  are  very 
close  to  global  optima  in  most  cases,  and  the  aver¬ 
age  percentage  excess  over  the  optimum  is  limited  to 
merely  within  1.12  percents  up  to  the  network  size  of 
20  nodes.  Moreover,  unlike  the  exact  method  by  inte¬ 
ger  programming  approach,  our  algorithms  can  handle 
very  large  network  sizes  and  hence  can  serve  as  the  ba¬ 
sis  for  benchmarking  other  developed  algorithms.  We 
believe  the  gained  insights  from  this  work  may  spur  the 


development  of  better  performing  localized  algorithms. 
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